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Abstract. We outline a discretization approach to determine the 
maximal number of mutually unbiased bases in dimension 6. We 
describe the basic ideas and introduce the most important defini- 
tions to tackle this famous open problem which has been open for 
the last 10 years. Some preliminary results are also listed. 
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1. Introduction 

This paper is based on the talk given by the second author at the 
International Conference on Design Theory and Applications, NUI, 
Galway, July 1-3, 2009. 

The notion of mutually unbiased bases (MUBs) constitutes a basic 
concept of Quantum Information Theory and plays an essential role 
in quantum-tomography [T5l [23] , quantum cryptography [UEEI2Q], the 
mean king problem [1] as well as in constructions of teleportation and 
dense coding schemes |22j. 

Recall that two orthonormal bases of C d , A = {ei, . . . , e^} and B = 
{fi, . . . , id} are said to be unbiased if, for every 1 < j, k < d, \ (ej, f^) | = 

—y=. A set Bo,...B m of orthonormal bases is said to be (pairwise) 

mutually unbiased if every two of them are unbiased. It is well-known 
( see e.g. [2j[5l[23]) that the number of mutually unbiased bases (MUBs) 
in C d cannot exceed d + 1. It is also known that d + 1 such bases can 
be constructed if the dimension d is a prime or a prime power (see e.g. 
dU 121 [El EEl HO [23]). Such a set of d + 1 MUBs in dimension d 
is called a complete set. If the dimension d = p^ 1 . . . p^ h is composite 
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then very little is known except for the fact that there are at least 
p" 3 + 1 mutually unbiased bases in C d where p^ 3 is the smallest of the 
prime-power divisors. Thus, the first case where the largest number of 
mutually unbiased bases is unknown is d = 6: 

Problem 1.1. 

What is the maximal number of pairwise mutually unbiased bases in 
C 6 ? 

Although this famous open problem has received considerable atten- 
tion over the past few years ([5j [7J [H [161 EH HI]), it remains wide open. 
Since 6 = 2 x 3, we know that there are at least 3 mutually unbiased 
bases in C 6 (see also [2U [16] for infinite families of MUB-triplets), but 
so far tentative numerical evidence [7J El [TOj [21] suggests that there are 
no more than 3, a fact apparently first conjectured by Zauner 



Conjecture 1.2. 

The maximal number of pairwise mutually unbiased bases in C 6 is 3. 

One reason for the slow progress is that mutually unbiased bases 
are naturally related to complex Hadamard matrices. Indeed, if the 
bases B ,...,B m are mutually unbiased we may identify each B[ = 
{ef \ . . . , e^} with the unitary matrix 



\h 



Hk,j 



<k,j<d 



i. e. the fc-th row of Hi consists of the coordinates of the fc-th vector 
of B\ in the basis B . (Throughout the paper the scalar product (., .) 
of C d is linear in the first variable and conjugate-linear in the second. 
Note also that for convenience of computer programming we use the 
unconventional definition that the rows of the matrices correspond to 
the vectors of the bases.) With this convention, Hq = I the identity 
matrix and all other matrices are unitary and have entries of modulus 
1/y/d. Therefore, the matrices VdHi have all entries of modulus 1 
and complex orthogonal rows (and columns). Such matrices are called 
complex Hadamard matrices. It is clear that the existence of a family of 
mutually unbiased bases Bq, ... , B m is thus equivalent to the existence 
of a family of complex Hadamard matrices y/dHi , . . . , yfd~H m such that 
for all 1 < j ^ k < m, VdHjH* k is again a complex Hadamard matrix. 
In such a case we will say that these complex Hadamard matrices are 
mutually unbiased. 

In particular, the existence of 4 MUBs Bq, 81,82, B3 in dimension 6 is 
equivalent to the existence of 3 mutually unbiased Hadamard matrices 
Hi, H 2 , H 3 . Therefore, in an attempt to prove that no collection of 4 
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MUBs exist in dimension 6 it is enough to prove that no triplet of 
mutually unbiased Hadamard matrices exist. This will be the core of 
our argument in this note. 

A complete classification of complex Hadamard matrices, however, is 
only available up to dimension 5 (see [H]) which allows for a complete 
classification of MUBs (see [9]). The classification in dimension 6 is 
still out of reach despite recent efforts [3J [191 EI]- This is one of the 
reasons for Problem II. II to be difficult. 

In this paper we outline a discretization approach that is likely to 
lead to the proof of Con] ect ure 1 1 . 2 1 in the near future. Once all the ideas 
are properly implemented in a computer code, an exhaustive search will 
be carried out to prove Conjecture 11.21 We will include all the basic 
definitions and ideas as well as some preliminary results. We remark 
here that a partial result of this approach has already been completed: 
in [TB] we assumed that the first Hadamard matrix Hi comes from 
the two-parameter Fourier family of complex Hadamard matrices, and 
we proved by discretization and an exhaustive computer search that 
in such MUB-quartet /, Hi, H2, H3 cannot exist. In this pa- 

per, however, we tackle the general case, so that we must consider Hi 
as any complex Hadamard matrix of dimension 6. This complicates 
matters quite considerably as the number of cases to check after the 
discretization increases by orders of magnitude. As an optimistic note 
let us recall here that the non-existence of a projective plane of order 
10 was also proved by an exhaustive computer search [TS] . 

2. Discretization 

The proof proceeds by contradiction, via a discretization scheme. 
Assume that there exists a collection of 4 MUB's in C 6 . Equivalently, 
there exist 6x6 complex Hadamard matrices A,B,C having all entries 
of modulus 1, such that the rows (and thus the columns) are complex 
orthogonal, and we have the unbiased condition: for any two rows 
u,v coming from different matrices we have \(u, v)\ = v6. (Recall 
that for the purposes of this note the rows of the matrices correspond 
to the vectors of the bases.) In such the orthonormal bases 

-jgA, -^B, -^C accompanied with the identity matrix Id correspond 
to a family of 4 MUB's. We assume that such matrices A,B,C exist 
and try to reach a contradiction. 

After multiplying rows and columns by appropriate scalars if neces- 
sary, we can assume that all coordinates of the first row and column of 
A are l's, and all coordinates of the first column of all other matrices are 
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l's (i.e. we assume that all vectors in the bases A, B, C have first coordi- 
nate 1, and the first vector in basis A is an all l's vector). All the other 
coordinates in the matrices are complex numbers of modulus 1, i.e. they 
are of the form e 2mp with p G [0, 1). We will use a discretization ap- 
proach. Let N be a positive integer, called the discretization parameter. 
We partition the interval [0, 1) into iV sub- intervals Iq, l[ N \ ■ ■ ■ , I^-i 
of equal length, i.e. J. = \j/N, (j + l)/N). (Other partitions are also 
possible, but this seems most convenient for programming.) Now, any 
entry e 27Tlp in any of the matrices A, B, C will be represented by the 
integer j if p G Ij N ^ (note that < j < N — 1). This means: whenever 
we see an entry j somewhere in a matrix then we conclude that the 
original phase p must lie somewhere in the interval /j . We have no 
more and no less information than this. We also agree that the first 
coordinate of each row will be represented by 0, keeping in mind that 
it represents exactly 1, without error (and not the interval Iq). 

In short: we will exclusively be dealing with row vectors of the form 
(1) u = (0, ji, j 2 ,j 3 , H,js) 

where < jk < N — 1 and the first coordinate represents 1 without 
error, while the other coordinates jk mean that the actual entry pk falls 
into the interval Ij^ ■ In notation, the original matrix will be denoted 
by A, while its representative integer matrix will be denoted by A. The 
entries of A will be denoted by p mt k, while those of A will be denoted 
by j m ,k- 

There are altogether N 5 vectors of the form ([T]). 

Also, there is a natural ordering among these vectors: u < v if 
and only if it is so in lexicographical order. We will use this ordering 
throughout this note. 

3. The search for the discretized Hadamard matrix A 

The matrix A is an integer matrix with first row and column con- 
sisting of 0's and the core of the matrix containing integers between 
and N — 1. We introduce the following definition: 

Definition 3.1. Given an integer matrix A with first row and column 
consisting of 0's and the core of the matrix containing integers be- 
tween and N — l, we will say that A is an iV-discretized representative 
of a complex Hadamard matrix if there exists a complex Hadamard ma- 
trix A with entries p m ^ such that p m ^ G Ij k . In notation A G HAD^, 
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where HADn denotes the set of N-discretized representatives of com- 
plex Hadamard matrices. 

The aim of this section is to describe an algorithm to efficiently search 
for all possible matrices A e HAD^- Upon strong numerical evidence 
PTj . it is conjectured that the manifold of 6 x 6 complex Hadamard 
matrices is 4-dimensional. Therefore we expect that the cardinality of 
HADjy will be approximately cN for some constant c. Nevertheless, 
the task of finding all possible A is daunting at first glance. There are 
N 25 possible iV-discretized matrices altogether, and we must select the 
ones belonging to HAD^- The number N 25 is of course astronomical 
even for N « 50, but we will see that with an intelligent approach the 
task can still be carried out. 

There are a few properties we can assume about A without loss of 
generality. We already assumed that the first row and column consist of 
0's. We can also assume that both the rows and columns are arranged 
so that they increase with respect to lexicographical order. This can be 
arranged by repeated permutation of rows and columns. (This is not 
entirely trivial because ordering the rows lexicographically can actually 
spoil such an ordering of the columns and vice versa. However, if one 
writes out the matrix entries row-after-row in one 36-long row vector, 
then it is clear that this vector will decrease lexicographically irrespec- 
tively of whether you make an ordering of the rows or the columns. 
Therefore such a repeated ordering of rows and columns will terminate 
in finite steps, and will produce a matrix such that both the rows and 
the columns increase in lexicographical order.) This automatically im- 
plies that the entries of the second row and second column are both 
monotonically increasing. This is a very convenient property, because it 
restricts the possibilities for the second row and column quite strongly. 

We can also assume that the second row is less than or equal to 
the second column in lexicographical order (this can be arranged by 
transposition of A if necessary). 

We must make use of the fact that the rows (and columns) of A 
are complex orthogonal to each other. The first row and column of A 
consist of 0's (representing the entry 1 in A, without error). Therefore, 
we have 5 unknown rows and columns of A. All of these rows and 
columns have the form ([I]). The orthogonality condition with the first 
row (and column) makes it natural to introduce the following definition: 

Definition 3.2. We will say that a vector u of the form (pQ) belongs to 
ORT N if there exist fc e I jk such that 1 + = 0. 
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Note that ORT^ is a "small" subset of all the vectors of form (pQ), 
containing only those vectors which represent vectors being orthogonal 
to the vector (1,1,1,1,1,1). Clearly, all rows and columns of A must 
belong to ORT^. Therefore it is very important to determine the set 
ORTn as precisely as we can. We achieve this by the following "check 
the descendants" method. 

Let u = (0, jx, j2, j3, j4, js), and let rj k denote the midpoint of the 
interval Ij k (the superscript N has been dropped from the notation for 
convenience). If u G ORTn then the trivial error bound (see Lemma 
3.1 in [16]) gives 




This is too crude, but we can iterate it to the "children" of u. Namely, 
assume that the numbers exist as in Definition [372J For each interval 
Ij k the value <pk must lie in either the left or the right half of Ij k . There 
are 32 choices, according to whether we consider the left or the right 
half of each interval Ij k . These choices are called the 32 "children" of 
it. Clearly, at least one of these children needs to satisfy ([2]) with |^ 
on the right hand side (and its own midpoints substituted to the left 
hand side, of course, instead of r Jfe ). If none of the children satisfy this, 
then u can be discarded. Of course we iterate this to grandchildren, 
and so on, down to 7-8 generations. The vector u survives this test if 
it has at least one surviving descendant in each generation. 

Remark 3.1. The set ORTn is clearly invariant under permutations 
of the last 5 coordinates ji, j 2 , J3, J4, J5. Therefore it makes sense to 
introduce the set ORT N rnon of vectors in ORT N with monotonically 
increasing coordinates. To save time, in the actual computer code we 
first find the vectors of ORT^ mon by the method above, and then we 
permute the last 5 coordinates to arrive at the set ORT^. 

Remark 3.2. There exists also an improved error bound (see Lemma 
3.2 in [IS]). It is somewhat slower to check by computer and it is 
reasonable to believe that we arrive at the same set ORT^ by applying 
either error bounds. 

Remark 3.3. We have implemented a computer code for selecting 
the set ORTn- For example, for N = 17 we have \ORT^\ = 58450, 
for N = 19, \ORT N \ = 82630, and for N = 53, \ORT N \ = 1875110. 
Experience shows that the set ORT^ is unexpectedly large if N is 
divisible by 2 or 3. Therefore, we have mainly restricted our attention 
to iV being a prime. 




(2) 




k=l 
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Remark 3.4. The optimal choice of N seems to be crucial for the 
success of the project. Clearly, if N is too small then the error bounds 
are not good enough and we will not reach a contradiction in the forth- 
coming steps (see Section H] below). However, if N is too large then 
the size of the sets ORT N and correspondingly HAD N will be far too 
large to be manageable. At present we believe that the optimal choice 
of N is around N ~ 50. 

Let us turn back to the construction of A. All rows and columns 
must come from ORT^, and they must be pairwise N — orthogonal in 
the following sense: 

Definition 3.3. We will say that the vectors u = (0, jx, 32, J3, ji, J5) and 

v = (0, mi, m 2 , m 3 , m 4 , m 5 ) are iV-orthogonal if there exist numbers <pk 
and ipk in the intervals Ij k and I mk , such that 1 + Ylk=i e 2m ^ k ~^ k ^ = 0. 

This property is clearly shift-invariant in the sense that it only de- 
pends on the values (ji - mi, . . . j 5 -m 5 ) modulo N. We can therefore 
take m\ = ■ ■ ■ = m 5 = and correspondingly v = (0, . . . , 0), (where 
the last 5 coordinates represent the interval Iq, of course) and define 
the set ORT epSt N as the set of vectors of the form (pP) which are iV- 
orthogonal to Vq. (The notation ORT epSt N indicates that the vector t> 
contains an "epsilon" of error, because the last 5 coordinates represent 
the interval Iq and not the exact number 1.) With this notation the 
shift-invariance means that u and v will be iV-orthogonal if and only 
if the vector (j 1 — mi, . . . J5 — m 5 )(mod N) is in ORT epsN . 

Having constructed the set ORT N previously, it is now easy to obtain 
ORT eps N . Indeed, by definition a vector u = (0, j\, . . . j 5 ) can only be 
iV-orthogonal to Vq if there exist numbers (f>k m the intervals Ij k and 
ip k in [0, ^), such that 1 + J2k=i e 2 ^ - ^ = 0. But then the numbers 
<Pk — ipk must fall in the intervals Ij k - ek where is either or 1, and 
hence the vector u e = (0,ji — ei, . . . , j 5 — e 5 ) is in ORT N . 

Therefore, ORT eps ^ will consist of all the vectors of the form u e = 
(0, ji +e u ...,j 5 + £ 5 ), where e k is or 1, and the vector (0, j u . . . , j 5 ) 
is in ORTn- 

Remark 3.5. Each u G ORT^ gives rise to 32 different u e above. 
One could therefore expect that the size of ORT epSj N will be nearly 32 
times the size of ORT^. This is not so, however, because there will 
be many coincidences. Experience shows that the size of ORT epSy jy is 
approximately 4 times the size of ORT N , regardless of the value of N. 
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Now we are ready to conduct a search for the possible matrices A. 
The first row and column are full of O's. We then build up the matrix 
with a row-by-column approach. We fit in the second row, then the 
second column, then the third row, then the third column, etc. At each 
step we must consider that: 

- each row and column must come from ORTjy. 

- each row (resp. column) must be lexicographically larger than 
any previous rows (resp. columns). In particular, the entries of the 
second row and column are monotonically increasing, i.e. they belong 

to ORTN )mon . 

- the second column must be lexicographically larger than or equal 
to the second row. 

- each row (resp. column) must be iV-orthogonal to any previous 
rows (resp. columns). This is equivalent to the fact that the pairwise 
differences of the rows (resp. columns) modulo iV must be contained 
in ORT eps N . 

- each row (resp. column) must be compatible with the already 
existing entries of the matrix (e.g. when we fit in the fourth row, then 
its first 3 coordinates are already fixed because the first three columns 
of the matrix have already been filled out previously). 

We have implemented a computer code which executes the search 
as described above. The running time is still reasonable, within 1-4 
days, depending on N. However, the number of selected matrices A is 
unexpectedly large. It is in the range of 10 9 — 5 • 10 10 as iV ranges from 
17 to 53. Let PREHAD^ denote the set of matrices obtained by this 
search. Clearly, HAD N C PREHAD N . 

Have we made all possible restrictions so as to list exclusively the 
matrices A belonging to HAD N ? In other words, is it true that 
HAD N = PREHAD N ? It turns out that this is not the case, and 
there is an important possibility for further pruning. Consider a ma- 
trix A G PREHADn. There are 25 non-trivial entries in A (the first 
row and column being trivial), all of which represent intervals Ij mk of 
length 1/N. Once again we can "check the descendants" of A. That is, 
we can take left or right halves of each 25 intervals L . , and therefore 

O Jm,k ' 

consider the 2 25 children of A. Obviously, at least one of these chil- 
dren need to satisfy stricter pairwise orthogonality conditions of rows 
and columns. If none of the children do, then A can be discarded, i.e. 
it does not belong to HADn- Of course, checking 2 25 children is very 
slow, but if one proceeds row-by-row then only a few thousand children 
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need to be actually checked. We have not rigorously implemented this 
step in our computer code. Nevertheless, preliminary results suggest 
that only a small fraction of the matrices in PREHAD^ will pass this 
test, i.e. HADn will be significantly smaller in size than PREHAD^. 
This is very important for the running time of the overall algorithm, 
as the size of HADjy should definitely be kept in the range 10 8 — 10 9 
even for N rs 50. 

4. Stage 2: vectors unbiased to A, and reaching a 

CONTRADICTION 

Let us fix a matrix A G HAD^. We want to prove that the pair 
(Id, A) cannot be extended by matrices B, C so as to meet all orthog- 
onality and unbiasedness conditions. The rows of B and C are of the 
form ([T]) and they must be "unbiased" to all six rows of A. Therefore, 
as a next step, we must obtain a list of all such vectors. 

Remark 4.1. We are actually free to use a different discretization 
parameter N' for the matrices B and C. It may well reduce the running 
time if we use optimal choices for N and N' . Experience shows (see 
[16] ) that it makes sense to choose N' considerably smaller than N. 
However, for the sake of simplicity we will keep N = N' throughout 
this note. 

As the first row of A is invariably (0, 0, 0, 0, 0, 0) (representing l's 
in the first row of A, without error) it makes sense to introduce the 
following definition: 

Definition 4.1. We will say that a vector u = (0, ji, j%, j'4, j$) be- 
longs to the set UB^ if there exist 0^ e Ij such that |l+X/fc=i e 2l7T ^ >k \ — 
VQ. We will say that u belongs to \JB^^ mon ij the coordinates of u are 
monotonically increasing. 

The set UB^ can be constructed in a similar way as ORT^. With 
rj k denoting the midpoint of the interval Ij k the trivial estimate gives 



This is too crude, of course, and the descendants of u need to be checked 
for some 7-8 generations. 

Remark 4.2. Once again, the set UBn is invariant under the permuta- 
tion of the last 5 coordinates ji, j'2, J3, ji, %■ Therefore, in practice, we 
first check monotonically increasing vectors only, and obtain UB^ tTnon . 
Then we permute the coordinates to obtain UB N . 





10 



P. JAMING, M. MATOLCSI, AND P. MORA 



Remark 4.3. The set UB^ is much larger than ORT^. This can be ex- 
pected because orthogonality of complex vectors induces two conditions 
(the real part and imaginary part both being zero) while unbiasedness 
only induces one condition. 

Remark 4.4. We have implemented a code for listing the set of vectors 
UBn- For example, for N = 17 we have \UB^\ = 479340, while for 
N = 19, \UB N \ = 764060. 

We will also need a set UB eps ^ which is analogous to ORT eps ^. 

Definition 4.2. We will say that the vectors u = (0, jx, j2, J3, J4, is) 

and v = (0, mi, m.2, m^, m 4 , m 5 ) are iV-unbiased if there exist numbers 
<p k andip k in the intervals Ij k and I mk , such that \ l+Ylk=i e 2 " 1 " 1 ^ - ^! = 
v/6. 

This property is again shift-invariant in the sense that it only depends 
on the values (ji — m 1; . . . j 5 — m 5 ) modulo N. We can therefore take 
rrii = ■ ■ ■ = m 5 = and correspondingly Vq = (0, . . . , 0), (where the 
last 5 coordinates represent the interval Iq, of course) and define the 
set UB epS) N as the set of vectors of the form (JI]) which are iV-unbiased 
to Vq. With this notation the shift-invariance means that u and v will 
be iV-unbiased if and only if the vector (ji — mi, . . .% — m§){mod N) 
is in UB epS:N . 

Let the vector v of the form ([1]) be any row of either B or C. As the 
first row of A is invariably (0,0,0,0,0,0) (representing the l's in the 
first row of A, without error), we conclude that v G UB N . Let a 2 , . . . a@ 
denote the last five rows of A. Then, by definition, the differences v — aj 
modulo ./V must belong to UB eps ^ for all j = 2, . . . 6. Let UB^ denote 
the set of vectors v which satisfy these conditions. The notation U B ^ 
reflects that these are the vectors which are "unbiased" to all rows of 
A. By what has been said above, all rows of B and C must belong to 
UB A . 

Remark 4.5. We have implemented a code to obtain the set UB A . 
Experience shows that the size of UB A is largely independent of the 
choice of A, and it has between 10 3 — 10 4 vectors as iV ranges from 17 
to 53. 

Having constructed UB A we must show that B and C cannot be built 
from these vectors satisfying all orthogonality and unbiased conditions. 

Consider the vectors in UB A and try to build the matrix B out of 
them. This means that we need to find 6 vectors b±, . . . bg such that the 
pairwise differences b k — b m modulo N all belong to ORT eps N . Counting 
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constraints and parameters one would expect that only a finite number 
of triplets of MUB's (Id, A, B) exists, and a MUB-pair (Id, A) can 
generically not be extended to a triplet (Id, A, B). This would give us 
hope that a contradiction is reached most of the times while trying to 
build B. However, recent results [lEJ IE] show that infinite families of 
MUB-triplets do exist. Numerical practice also shows that the matrix 
B can indeed be built from the vectors of UB^ for all A. Therefore, we 
do not get an immediate contradiction. Instead, for each B we must 
go on and select the vectors UB^ ^ which are unbiased to all rows 

of A and B, and we must try to build a matrix C out of the vectors 
UB^ g. The contradiction is reached only at this point. That is, if N is 

large enough the matrix C cannot be constructed from UB^ g to meet 
all orthogonality conditions. Experience shows that N must be larger 
than 30 to reach a contradiction. This part of the project is currently 
under implementation. It would be desirable to reach a contradiction 
for each A within a few seconds of computing time. 

For the overall success of the project two tasks need to be considered 
in the near future. One is the implementation of the ideas of the last 
paragraph of Section [3]to bring down the number of possible A's to the 
region 10 8 — 10 9 . The other is to reach a contradiction for each A within 
a few seconds of computing time by not being able to construct B and 
C . A nice feature of the overall project is that once the algorithm is 
completed, it is very easy to distribute the calculations among several 
hundreds of computers, and thus reducing the running time by 2-3 
orders of magnitude. 

Finally, we remark that the entire discretization procedure described 
above has already been completed in [16] in the restricted setting 
when A is assumed to belong to the Fourier family F(a, b) of complex 
Hadamard matrices. 

Theorem 4.6. [Theorem 1.4 in |16j ] 

None of the pairs (id, F(a, b)) of mutually unbiased orthonormal bases 
can be extended to a quartet (id, F(a,b), B,C) of mutually unbiased 
orthonormal bases. 

In that case we used the discretization parameters N = 180 for 
A and N' = 19 for B and C. Due to some well-known equivalence 
relations only a few hundred possible discretized matrices A needed to 
be considered, and a contradiction was quickly reached for all of them. 
The documentation of that search is available at [25] . The difficulty in 
the general case is that the number of matrices A becomes very large 
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if iV is chosen large, while if N is small then a contradiction is reached 
very slowly (or not reached at all!) in the second stage of the search. 
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